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ABSTRACT 

We have undertaken a project to investigate the host galaxies and environ- 
ments of a sample of quasars at 2; ~ 4. In this paper, we describe deep near- 
infrared imaging of 34 targets using the Magellan I and Gemini North telescopes. 
We discuss in detail special challenges of distortion and nonlinearity that must 
be addressed when performing PSF subtraction with data from these telescopes 
and their IR cameras, especially in very good seeing. We derive black hole masses 
from emission-line spectroscopy, and we calculate accretion rates from our Ks- 
band photometry, which directly samples the rest-frame B for these objects. We 
introduce a new isophotal diameter technique for estimating host galaxy lumi- 
nosities. We report the detection of four host galaxies on our deepest, sharpest 
images, and present upper limits for the others. We find that if host galaxies 
passively evolve such that they brighten by 2 magnitudes or more in the rest- 
frame B band between the present and z=4, then high-z hosts are less massive 
at a given black hole mass than are their low-z counterparts. We argue that 
the most massive hosts plateau at < lOL*. We estimate the importance of se- 
lection effects on this survey and the subsequent limitations of our conclusions. 
These results are in broad agreement with recent semi-analytical models for the 
formation of luminous quasars and their host spheroids by mergers of gas-rich 
galaxies, with significant dissipation, and self-regulation of black hole growth and 
star-formation by the burst of merger-induced quasar activity. 



Subject headings: galaxies: evolution, high-redshift — quasars: general 
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Introduction 



In the past decade, we have begun to understand the important role that black holes 
play in galaxy evolution. Observations suggest that supermassive nuclear black holes are 
likely present in nearly all normal galaxies, and that black hole mass i s cor related with 



host galaxy bulge mass (Kormendv fc Richstonel IQQsl: Magorrian et al. IQQsh and stellar 
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These results have implications for the evolution of luminous, high- 



redshift quasars. If a galaxy already had a supermassive black hole early on, then according 
to the local black- hole/bulge relation, it must by today be one of today's most massive 
galaxies. 

In the context of the ACDM framework for hierarchical structure growth, specific predic- 
tions can be made for the evolution of quasar hosts an d their environments through cosmic 
time (IMo fc Whitdl2002l : iKauffmann fc Haehneltll2002l ). If dissipationless gravitational col- 
lapse of cold dark matter were the only process at work, then one would expect the ratio 
of black hole mass to stellar spheroid mass, Mbh/M^, to be roughly constant as spheroids 
merge and their nuclei coalesce. However, luminous quasars like the ones in current samples 
of quasars at z > 4 are likely the product of ma. i or mergers of gas-rich disk galaxies of com- 



parable mass (|Hopkins et al.l 12005 



Crotonll2006l : IPi Matteo et al.ll2008l : iHopkins et al.ll2008 



Somerville et al.ll2008l ). The central black holes merge, and merger-induced gas accretion re- 
sults in a burst of quasar activity. Quasar radiative energy and winds eventually halt further 
mass accretion and clear out the cold gas, halting star-formation. The merger disrupts the 



^Based on data obtained with the 6.5 meter Baade Telescope of the Magellan Telescopes, located at the 
Las Campanas Observatory, Chile. 

^Based on observations obtained at the Gemini Observatory, which is operated by the Association of 
Universities for Research in Astronomy, Inc., under a cooperative agreement with the NSF on behalf of the 
Gemini partnership: the National Science Foundation (United States), the Particle Physics and Astronomy 
Research Council (United Kingdom), the National Research Council (Canada), CONICYT (Chile), the 
Australian Research Council (Austraha), CNPq (Brazil) and CONICET (Argentina). 

■^Based in part on data taken at the W. M. Keck Observatory, which is operated as a scientific partnership 
among the California Institute of Technology, the University of California, and NASA, and was made possible 
by the generous financial support of the W. M. Keck Foundation. 
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original gaseous disks, and the result is a low angular-momentum spheroid of stars that subse- 
quently evolves passively. Semi-analytical models for these processes, along with numerical 
simulations of the collapse of cold-dark-matter halos, are successful in reproducing many 
observations of galaxies, galaxy clusters and quasars. In these mo dels, high-z quasars are 



expec ted to have less luminous hosts than their low-z counterparts (IKauffmann fc Haehnelt 



2OOOI ). w ith the ratio of black hole mass to st ellar spheroid mass, Mbh/M^,, decreasing with 
redshift (1Croton|[2006l : ISomerville et al.lboosh . 



By necessity, models for quasar host evolution rely on semi-empirical prescriptions for 
key physical processes, because the resolution of numerical simulations cannot follow all of 
the crucial physics from the spatial scales of galaxy clusters down to galactic then to atomic 
scales. Direct observations of high redshift quasar hosts such as the study described here 
provide one interesting empirical check on the overall validity of the theoretical picture of 
hierarchical galaxy and black hole formation and evolution. 

Detecting the host galaxy "fuzz" is technically challenging at high redshift however 
because it appears small and faint compared to scattered light from the nucleus in the wings 
of the point spread function (PSF). Ideally, we would study the fuzz in the rest-frame near-IR, 
which would both highlight the mass-tracing st ellar populations of the hosts and provide the 
best possible galaxy-to-nuclear light contrast (IMcLeod &: Riekd Il995l ). For high-z objects, 
that would mean observing in the mid-IR, but there are not yet telescopes with the necessary 
combination of sensitivity and angular resolution to make such observations feasible. Most 
high-z host studies so far have therefore used near-IR imaging. 

At z ~ 2 — 3, the first near-IR imaging of handfuls of objects using 4m-class tele- 



scopes produced detect i ons only in the case of radio-loud (RL) quasars (ILehnert et al.l Il992 



Lowenthal et al. I Il995l : lAretxaga. Terlevich. fc Boyld Il998l : ICarballo et al.l Il998l ) . Fuzz was 
subsequent ly seen around a few radio-quiet (RQ) quasars using adaptive optics (A O) on 4m 
telescopes ( Aretxaga et al. 1998 : Hutchings et al. 1999 : Kuhlbrodt et al. \ 2005 ). and AO 
on the Gemini North 8m yielded only one of 9 hosts aX z ^ 2 (jCroom et al.ll2004l ). The well- 
characterized and stable PSF of NI CMOS allowed successful detections of larger samples 



of RQ hosts in this redshift range (iKukula et al. I l200ll : iRidgway et al. 1 12001 
2006[ ). which led to the first tantalizing comparisons with hierarchical models 



Peng et al. 



Peng et al. 



(I2OO6I ) suggest that at 2; = 2, Mbh/M^, was several times larger than it is today, and that 
hosts were n ot yet fully for r ned, th ough those results have been the subject of some debate. 
For example iFalomo et al.l (120081 ) have recently used AO imaging on the VLT to detect 
the hosts of three luminous quasars at 2 < 2; < 3, which they use to help argue for pas- 
sively evolving ~ 5L* elliptical hosts. Moreover, Malr nquist bias in high-z samples would 
skew detections towards bright quasars in small hosts (ILauer et al. 1120071 : iTreu et al.l 120071 : 
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Woo et al.ll2008l ): we discuss this effect further in §7.61 



To provide more leverage for probing the hierarchical models, we would like to measure 
host properties at 2; ~ 4. So far, host de tections have been claimed for just a few objects in 
this range, all with unknown radio type. iPeng et al. I ( l2006l ) used NICMOS to measure host 
magnitudes and sizes for two gravitationally lensed quasars at z = 4 . 1 and 4.5. Hutchings 
has used Gemini to observe seven targets with 2; ~ 5 (lHutchingsll2003l . 120051 ): unfortunately, 
the results are hard to interpret in light of the distortion and nonlinearity effects that we 
have found to be significant for the same instrument (see below). 

With the improvement in sensitivity provided by 6-8m class telescopes, we have begun 
to expand the sample of quasars imaged at high z. Our long-term program includes new 
multi-wavelength imaging to search for hosts and to characterize environments; new and 
archival visible spectroscopy to be used for virial mass estimates of the central engines from 
emission line properties; and modeling of the accretion disks and hot coronae using data 
from IR to X-ray, including new X-ray data from ongoing Chandra and XMM observing 
programs. 

In this paper, we describe a sample of 34 z ~ 4 quasars that we have imaged in the near- 
IR. We present the observations and report on the search for hosts. The environments will 
be described in a subsequent paper (Bechtold & McLeod 2010, in preparation). We adopt 
the cosmology = 0.7, Q^a = 0.3, Hq = 70 km s~^ Mpc~^ throughout. All magnitudes are 
reported in the Vega system, unless otherwise noted. 



2. The Sample 



We have observed a sample of quasars selected to have redshifts in the range 3.6 ^ z < 
4.2. The sample is listed in Tabled] with names as given in the NASA/IPAC Extragalactic 
Database (NED). The redshift range was chosen so that the 4000A break falls between 
observed H- and i^'-bands, so that broad-band colors give maximum leverage for estimating 
photometric redshifts and stellar populations. Out of the ~ 300 quasars known in this 
interval when we began the project, we observed a randomly-chosen sub-sample of 34 objects, 
yielding a median < z >= 3.9 and spanning a range of magnitude as shown in Fig. [H We 
plot our sample against the ~ 1600 quasars at these redshifts listed i n the most recent Sloan 
Digital Sky Survey (SDS^I) Quasar Catalog (jSchneider et al. 1120071 ). 



''Funding for the Sloan Digital Sky Survey (SDSS) has been provided by the Alfred P. Sloan Foundation, 
the Participating Institutions, the National Aeronautics and Space Administration, the National Science 
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We are observing first in the i^'-band, which samples the rest frame B for these objects. 
The median observed total nuclear Kg = 17.2 for the quasars in our sample corresponds to 
Mb = —26.9 (Vega magnitudes), similar to the local luminous quasar 3C273. 



2.1. Radio Loudness 



To characteri ze the radio prope rties of our sample, we adopt the definition of radio 
loudness given in llvezic et al.l (120021 ). wh o analyzed the rad io properties of SDSS qusars. 
They define the apparent AB magnitude f Oke &: Gunn 19831 ) at 1.4 GHz as 



-2.5logio 



int 



3631J?/ 



where Fint is the integrated 20cm radio flux measured from a two-dimensional Gaussian fit 
to the radio source. The radio-to-optical flux density is then defined as 



0.4(^^^ - 1) 



where i^^ i s the AB magnitude at Sloan i in the continuum. With these definitions 

1—4 and radio-quiet quasars 



Ivezic et al.l ( 120021 ) find that radio-loud quasars have Ri = 
have Ri <1. 

The radio properties for the quasars in our sample are given in Table [2l We compiled 
the optical magnitudes from a number of sources. When available, we adopted the values 
for i"^^ published by members of the Sloan consortium; references are listed in Table [21 For 
most other objects, we used the 2^-^-band photometry given on the SDSS web site (DR6). 
For two objects we measured photometry from archival HST images. For these and a few 
other cases, the only available photometry was from the literature in other filters, which we 
transformed to i^^ using the zero-points given in the NICMOS web site unit converter, and 
assuming a quasar spectral energy distribution in the form of a power law 

(X v"" where the spectral index a = -0.44 for 1216A < Xrest < 5000A 



which is the aver a ge de rived from the SDSS quasars by IVandenBerk et al. I ( l200ll ) (see also 
Pentericci et al. I J2003h l 



Foundation, the U.S. Department of Energy, the Japanese Monbukagakusho, and the Max Planck Society. 
The SDSS Web site is |http: / /www.sdss.org/ 
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We corrected t he op tical magnitudes for Galactic reddening, using the E{B — V) given 



bv lSchlegel et al. I fll998f ) as tabulated in NED. The transformat ions A,. = E{B - F)/2.751 



and Ai = E{B - l^)/2.086 were adopted fISchneider et al. Il2003[ ). 



For most of the sample quasars, the mos t sensitive radio da ta comes from the Faint Im- 
ages of the Radio Sky at Twenty-cm Survey (IBecker et al.lll995l . FIRST), which we accessed 
through NED. If the quasar did n ot fall in one of the FIRST fields, we used the data from 
NRAO VLA Sky Survey (NVSS; (ICondon et al.lll998l )). In a few cases the most sensitive 
radio data had been reported in targeted searches in the literature. For objects detected in 
FIRST, we adopted the FIRST catalog integrated flux. For others, we derived a 2a upper 
limit from the root-mean-square fluxes in the maps downloaded from NED. 

Of the 34 sample quasars, 16 are radio-quiet, 5 are radio- loud, 4 have no radio data, 
and 9 have radio data which are not deep enough to know whether or not the quasar is 
radio- loud or radio-quiet. 

Going into our survey, we wanted to test whether radio-loud quasars have different host 
galaxy properties than the radio-quiet majority. We therefore tended to give priority to 
observing quasars which we knew to be radio-loud, since they are rare and we knew that it 
would be difficult to get a statistically large sample of them. In the end, at least 5 quasars in 
the sample of 34 are radio-loud, compared to approximately 1 expected to be radio-loud, had 
we observed a sampl e representative o f the radio properties of the bright quasar population 
at z ~ 4 as a whole (I Jiang et al.l 120071 ). 



2.2. Black Hole Mass Estimate 

We estimated the black hole mass Mbh for the quasars in our study from emission- 
line spectroscopy. As described below, we measured the full-width-half-maximum (FWHM) 
value of the broad CIV emission line and the quasar UV continuum luminosity from new and 
existing spectra. We then used these to compute black hole mass according to the relation 



log 



10 
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Mo 



log 



10 



FWHMciv y f XL 



^A1350 



1000km s-^ J yW'^erg s-^ 



0.53 



+ 6.66 



from lVestergaard fc Peterson! (l2006l ). who find that the UV continuum luminosity Lai45o can 
be freely substituted for Lai35o- From this we derived 
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Here, AB 



1450 



—2.5logiofiy — 48.6, where /j, is the (reddening-corrected) continu um flux in 
measured at an observed wavelength of A = 1450A(1 + z) as in lFan et al. 
( I2OOII ). DM is the (luminosity) distance modulus. 



The AB1450 magnitudes were compiled from the literature or measured by us as shown 
in Table [31 Where available, we adopt the reddening-corrected AB1450 tabulated by members 
of the Sloan consortium, who give values based on spectrophotometry of the quasars at a 
rest frame wavelength of 1450A. For 10 of the objects, we measured the continuum fluxes 
ourselves from spectra obtained from the SDSS Skyserver. For the objects for which no flux- 
calibrated spectra were available, we used the i^^ magnitudes from Table [2] and transformed 
to AB1450 assuming a = —0.44 as described in §2.11 For objects with large CIV equivalent 
widths that contaminate the broadband measurements, the AB1450 derived from photometry 
will be systematically bright. A comparison of the spectroscopically-derived AB1450 to the 
photometrically-derived one for the SDSS objects shows the former to be fainter on average 
by 0.4 ± 0.3 mag. 

For most objects, we measured the FWHM of the C IV emission lines given in Table 
E] using spectra from the SDSS Skyserver or electronic versions of published spectra from 
several authors who kindly made them available. In a few cases, we digitized published 
spectra using "Plot Digitizer" software. We also carried out new long-slit optical spectroscopy 
for flve targets in the sample, including one for which no other spectroscopy is published, 
[VII95]2125-4529. For the new obse rvations, we used the DEIMOS spectrograph on the 
Keck-II telescope JDavis et all l2003h on the nights of 2008 Oct 24 and 2004 Oct 12 and 
13. The objects were observed through a 0.7 arcsec wide slit with the 1200 1/mm flrst 
order grating, resulting in a dispersion of 0.33 A/pixel. A GG495 fllter was used to block 
second order light. Exposures were 600-900 seconds, mostly through clouds or at twilight. 
We reduced the data using the DEEP2 project IDL reduction pipeline, which flatflelded, 
sky-subtracted, wavelength-calibrated and extracted the spectra as described in the DEEP2 



webpage, http: / / astro. berkeley.edu/^cooper/deep/spec2d/| Spectra are shown in Figure [2l 



To derive the CIV line width, we subtracted a local continuum flt, derived by fltting a 
linear curve through the spectra in rest wavelengths 1425A to 1500A and 1760 A to 1860 
A. We replaced absorption features with an interpolated continuum estimate, and then flt a 
gaussian to the C IV emission line. 



The quality of the C IV line proflles for the quasars in our sample ranged from very high- 
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signal-to-noise examples with easy-to-define continua, to barely detected lines in discovery- 
quality spectra. Moreover, the redshifts of our targets shift the C IV line to wavelengths with 
strong telluric absorption and night-sky emission features that are difficult to calibrate out 
completely. Some lines probably are suppressed by undetected absorption features intrinsic 
to the quasar. As many authors have noted, quasar emission lines are non-Gaussian, in the 
sense that they have "pointy" peaks. Some C IV lines in our sample were also significantly 
asymmetric. 

For these reasons, we measured the FWHM values by hand, using IRAF'^ splot. In 
cases where part of the line profile was very noisy, we measured the half-width of the better 
side of the profile and doubled it. In a few of the spectra with very good signal-to- noise, 
there is clearly a narrow (2000-3000 km FWHM) component, and broader wings (10,000- 
15,000 km s~^ FWHM). For several objects, we estimated two values of the line width, one 
for each component; both are listed in Table [31 For [VCV96]Q2133-4625, the C IV profile is 
so noisy, possibly because of an absorption trough, that it was impossible to derive a FWHM 
from the published spectrum. 

Our resulting black hole mass estimates are given in Ta ble Bl The systematic uncertain 
ties in these estimates for black hole mass are well-known ( Wandel et al. Il999l: ICoUin et al. 



2002 



2007 



Dietrich fc Hamann 



Kellv fc Bechtold 



2004; 



2007 



Kaspi et al 



McGill et al. 



20051 : IVestergaard fc Petersonll200( 



Netzer et al. 



20081 ). The primary assumption is that the CIV 



emitting gas is in virial equilibrium with the central black hole mass, and is located at a 
radius that scales with luminosity. For the luminous quasars in our sample, this means 
extrapolating fro m the relations teste d in emission-line regions studied with reverberation 
mapping locally (jPeterson et al.ll2004j ). Further, the bolometric luminosity of each quasar 
is assumed to be a constant multipl e of the A1450 continuum luminosity, which certainly is 
not the case (e.g. iKellv et all tooi )). 

In Figure [3] we plot the black hole masses of the quasars in our samp le along with thos e 
for the ~ 1600 SDSS quasars in this redshift range recently tabulated by lShen et al.l (120081 ). 
For the 6 objects in common, our black hole mass estimates generally agree within ^ 0.2dex. 



'"^IRAF (Image Reduction and Analysis Facility) is distributed by the National Optical Astronomy Ob- 
servatories, which are operated by the Association of Universities for Research in Astronomy, Inc., under 
contract with the National Science Foundation. 
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2.3. Accretion Rates 



We combined tlie black tiole masses witli ttie K-band observations described below 
to calculate the quasar mass accretion rates. Because the observed K-band samples the 
rest-frame B-band, the K-band magnitude allows us to compute a B-band luminosity inde- 
pendently of spectral shape. This avoids the errors that result when one must extrapolate 
from optical photometry to the rest-frame B assuming a sp ectral index a. We apply a B- 
band bolometric correction factor of 10.7 (lElvis et al.lll994l ) and we compare the resulting 
bolometric luminosity to the Eddington luminosity computed from the black hole mass via 
LEdd = 3.3 X 1O^(Mbh/^0) Lq. We have assumed that all of the rest-frame B-band light 
can be attributed to the nucleus, which is a reasonable estimation for such luminous objects. 
The resulting accretion rates as fractions of Eddington, Lboi/LEdd, are tabulated in Table [31 

The median L^oi/LEdd for the sample is 0.47 ±1.6 (Icr), and the minimum value is 
0.1. These rates ar e good matches to those inferred from studies of host galaxies locally; 
McLeod fc McLeodl (I20011) found that the most luminous local quasars radiate at > O.IL^;;^;^, 
and iFloyd et al. I (120041 ) deduce a median rate 0.47 for the most luminous local quasars. 



The calculation of accretion rates yields a ha ndful of quasars with super-Eddington rates. 
Of these, BR2212-1626 is gravitationally lensed (jWarren et al.ll200ll ) and so the continuum 
luminosity, which we have not corrected for gravitational magnification, is overestimated. 
Since both Mbh and L/LEdd oc L°-^, both quantities are also overestimated. For BRJ0529- 
3553, we have only a discovery quality spectrum, and the C IV line width is very uncertain. 
For 5 others which have L/LEdd >> 1, the C IV profile has good enough signal-to- noise to 
detect a distinct narrow and broad component. If the FWHM of the broad component is 
used, very large black hole masses, and sub-Eddington accretion rates are implied. Detailed 
modeling of the quasar spectra energy distribution and higher quality spectra of all targets 
would improve the estimates of black hole mass and accretion rate. We do not list the 
statistical errors for Mbh and accretion rate in Table [31 because these numbers are dominated 
by systematic uncertainties and the simplifying assumptions described above. 

Excluding the L/LEdd > 1 objects , the m edian rate becomes 0.41 ±0.3, consistent with 
the distribution plotted by lShen et al.l (120081 ) for z > 3 SDSS quasars. 



As a second way to estimate black hole masses for our sample quasars, we assume that 
all of the quasars are radiating at OALEdd, with bolometric luminosities determined from 
Mb as above. The derived values are then the minimum plausible black hole mass that the 
nuclei could have to be emitting at the luminosity observed. These values are listed in Table 

El 
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3. Near-IR Imaging Observations 



We have obtained deep, near-IR images of 34 quasars over tlie period 2002 September 
- 2005 January at tlie Magellan I 6.5m and Gemini North 8m telescopes. We have observed 
each field in Kg, with 5 a lso observed in -H^ or H^. Most of the objects (26) were observed 
with Magellan's PANIC jMartini et al. Ibood ). a 1024 x 1024 HgCdTe array with a pixel 
scale of 0'.'125 and a field-of-view (FOV) o f 128^^ Before PANI C was installed, we imaged a 
few objects (6) with the old ClassicCam (jPersson et al. Ill992l ). a 256 x 256 HgCdTe array 
camera yielding a FOV of only ~ 30^^ per exposure. The rest of our targets (7) were observed 
on Gemini with NIRI Jflodapp. K. W. et al. Il2003h . a 1024 x 1024 InSb array operated at 
f/6, yielding a pixel scale of 0'.'116 and a FOV of 119" per exposure. The observations are 
summarized in Table [TJ 

With all three instruments, we observed using a 9- or 25- point dither pattern of short 
exposures (10 — 30 sec each, repeated 1 — 3 times per dither position). The times were chosen 
to keep the quasar images in the nominal linear range, with repeats limited to ensure fair 
sampling of sky variation. The dithers were typically repeated for half a night, yielding up 
to a thousand frames per field and average total on-source times of 3 hours for PANIC and 
NIRI. With NIRI and PANIC, the FOV was big enough for us to use stars on the quasar 
frames to measure the PSF. Due to ClassicCam's smaller FOV, we had to alternate quasar 
dithers with dithers on a nearby star to sample the PSF. The average on-source time for 
ClassicCam was thus shorter, only about 2 hours, resulting in shallower images. As we 
describe in the sections below, the ClassicCam images turned out to be of very limited use 
for the host searches. However, we include them in this paper both because they remain 
somewhat useful for investigations of the near environments of the quasars, and to illustrate 
the difficulties of using out-of-field stars for PSF subtraction. 

We were fortunate to have excellent seeing for many of the observations, with the final, 
combined Kg quasar images from PANIC and NIRI having full-width half-max (FWHM) 
ranging from 0'.'32 — 0'.'66 with median value < FWHM >= 0'.'43. ClassicCam images were 
worse. Photometric calibration was done for the PANIC and NIRI images using the 2MASS 
stars (usually 1—3) found on each combined quasar frame. We also cross- checked these values 
against observations of Persson faint IR standards (jPersson et al. Ill998l ). and estimate that 
the photometric calibration is accurate to ~ O.lmag. For the ClassicCam images, whose 
FOV are too small to contain 2MASS stars, we used only the Persson et al. standards. 



- 11 - 



4. Data Reduction 

For all three instruments, we reduced the data using standard techniques in IRAF 
with its add-on packages gemini/niri and panic, the latter kindly provided by Paul Martini. 
The data from each instrument were handled somewhat differently, with flats made from 
twilight exposures, fiat lamp exposures, and object frames for PANIC, NIRI, and ClassicCam 
respectively. Sky frames were generally made by median-filtering 9 — 10 dither positions after 
masking out sources. For the NIRI images, persistence proved to be a significant problem, 
which we solved to our satisfaction by including in each frame's bad pixel mask the object 
masks from the 2 previous frames. 

The PANIC and NIRI detectors are known to be nonlinear by ~ 1% at 15,000ADU 
and 11,000ADU respectively, and we kept our quasar+sky counts well below these limits. 
Still, for the tight tolerances in this project we needed to take some care with the linearity 
correction. For PANIC we used fiat lamp exposures of varying lengths to determine our 
own second-order correction which differed from the nominal pipeline correction by 0.5% at 
15,000ADU. The NIRI pipeline does not offer any nonlinearity correction and we lacked the 
data to determine our own. We discuss the residual effects of nonlinearity in ^ 

In the course of our analysis we detected a geometric distortion in the PANIC images. 
The distortion was visible as a radial stretch in contour plots of stars taken around an 
image. Paul Martini gave us a second-order geometric distortion correction derived from the 
PANIC optical prescription, which we then implemented in the pipeline. We note that the 
NIRI pipeline does not offer a distortion correction, though distortion proved to be an issue 
there as well. We discuss the implications further in §3 

The hundreds of reduced frames for each quasar were magnified by a factor of two, 
aligned on the quasar centroid, and combined after rejecting a handful of frames deemed 
bad because of bias level jumps in one quadrant or poor flattening. The deepest NIRI Kg 
images reach a surface brightness limit of Kg = 22.9mag arcsec"^ (measured as a la pixel- 
to-pixel variation). The median value is 21.7 mag arcsec"^, typical for the PANIC images, 
while the ClassicCam images are more shallow. A typical PANIC image is shown in Fig. 
m where the FOV corresponds to ~ 1.3 Mpc at z = 4, well-suited for studies of quasar 
environments (see Bechtold & McLeod 2009). 

5. PSF Characterization 

Any search for host galaxies is only as good as the characterization and removal of the 
nuclear point source. We followed the traditional practice of selecting "PSF stars" from each 
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image, and using them in model fits to the quasar images. However, in the course of our 
analysis we discovered some subtle effects of residual distortion and nonlinearity. Because 
we have not seen these issues addressed in other high-z host searches, including ones also 
done with NIRI, we discuss them in some detail here. For a recent look at th e PSF perils 
that host galaxy studies might encounter even with HST, see iKim et al.l (120081 ). 



5.1. Geometric distortion, or What to Do When The Seeing is Too Good 

Beginning with PANIC, our experiments with multiple PSF stars showed that poorer 
fits tend to result when the PSF star is farther from the quasar. Even though we had 
performed a distortion correction on the PANIC images, we found that a small residual 
geometric distortion compromised the fits. The distortion we detected would be insignificant 
(and indeed not noticeable) for most projects, with camera optics generally designed to create 
instrumental PSFs small compared to the seeing size. In our case, however, the excellent 
seeing and tight tolerances required for high-z host detection made the distortion apparent. 

To improve the fits, we were able to effect a higher-order distortion correction by re- 
centering each quasar's hundreds of frames on PSF stars to create "PSF frames" for each 
target, as suggested to us by Brian McLeod. With distorted images, the pixel scale at the 
edges is different than that near the center. Therefore, when shifting and combining frames 
from different dither positions, the shifts computed based on the objects near the center (in 
this case quasars) will be the wrong number of pixels to align the objects near the edges, 
yielding a combined image with stretched edges. The recentering technique to help correct 
for this works as follows. First, as described above, we align the hundreds of images to the 
quasar centroids, and combine them to create a "quasar frame." From this we extract a 
postage-stamp image of the quasar to use for fitting. We then start again with the same 
hundreds of images and align them this time to the centroid of a particular PSF star, and use 
these new shifts to combine them to make the "PSF frame" for that star. We then extracted 
a postage-stamp image of the PSF star from the latter frame for use with the fitting. An 
example is shown in Fig. [51 We used this technique with good success on most of the PANIC 
images. Examples of fits performed with and without our re-centering technique are shown 
in Fig. [61 

The NIRI images also suffer from distortion that proved significant for this project. No 
distortion correction is used in the NIRI pipeline. We applied our re-centering technique 
and found that it did improve the NIRI fits considerably, but PSF-PSF tests (where we 
subtracted PSF stars from each other) showed that residual distortion remains in the K image 
of SDSSJ012019.99+000735.5 and the H image of BRI0241-0146. For these two images, the 
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only PSF star is far from the quasar. 

The ClassicCam images provided their own set of challenges because the PSF stars were 
observed alternately with the quasars, and inevitable seeing variations resulted. We were 
able to obtain more reasonable results in most cases by rejecting frames according to the 
seeing, so that the resulting quasar and PSF star frames had the same FWHM. However, 
the ClassicCam results are never as satisfactory or robust as the Panic and NIRI results, 
and will be more useful for studies of the quasar's near environment than for host detection. 

5.2. Nonlinear ity 

Unfortunately, we also discovered that the near-IR images exhibit a small nonlinearity 
even after the nominal correction has been applied. This effect was more subtle than the 
distortion and became apparent only under scrutiny of the ensemble of data for our many 
objects. Such an effect could easily have escaped our detection in a study with fewer objects. 
We noticed that our best fits were found for PSF stars with similar brightness to the quasar. 
PSF stars brighter or fainter than the quasar could leave compact central emission or, more 
insidiously, rings that mimicked host galaxies in the difference images. Star-minus-star 
experiments performed on multiple PSFs from the same image confirmed our suspicion. 

This is difficult to illustrate with observed stars because of the residual distortion dis- 
cussed above. However, we investigated this further by simulating observations of stars of 
different magnitudes and fitting and subtracting them after applying different plausible lin- 
earity corrections. For example, we have used the nonlinearity curve for PANIC shown in 
Fig. El to generate the suite of stars shown in Fig. [Hlwith radial intensity profiles given in Fig. 
[9l The counts for the fainter stars were chosen to keep the detector within the range where 
the response is approximately linear, as is done for the observed targets. For the brightest 
stars, we allowed the brightness to enter the nonlinear (but not close to saturated) regime. 
For the PANIC response, this corresponds to ~15000 counts. At this level, the difference 
between the plausible prescriptions for the linearity correction amounts to < 0.5%. 

When we generated stars using one prescription and then "corrected" them for nonlin- 
earity using another, the 2D residuals in star-star tests were clearly positive. In other words, 
the uncertainty in the linearity correction can lead to spurious detections when scaling and 
subtracting point sources that differ in flux. However, in the cases we tried, the spurious 
residuals were distinguished either by unphysically compact sizes (FWHM less than the im- 
age FWHM) as shown in Fig. [9l or else by donuts that could be mistaken for over-subtracted 
hosts. The latter did not extend past a diameter of = 2.5 FWHM, as illustrated in Fig. 
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m 



5.3. Implications 

Our results call into question the traditional approach of selecting "PS F stars. .. chosen t o 
be as bright as possible without encountering detector saturati on effects" (Hutchingsl 2005 ). 



We have used this approach ourselves for low-z quasars (e.g. iMcLeod fc Riekel JlQQsh ) 



so 

that noise in the PSF wings is scaled down during the fitting process. However, our current 
analysis suggests that for PANIC and NIRI at least, a more robust practice is to chose PSF 
stars whose brightness is similar to the quasar, and whose positions are as close as possible. 
Which criterion takes priority might depend on the instrument and the observing conditions. 

Distortion corrections and linearity corrections are essential, but not sufficient. The PSF 
star re-centering technique described above provides a higher-order distortion correction, but 
a possible added complication is that the accuracy of the registration can be dependent on the 
brightness of the stars. In addition, the characteristics of spurious residuals are dependent 
on the weighting process used during normalization of the PSF (e.g. normalize to the flux 
in the central few pixels, use the whole source for the fit, weight the fit by flux, etc.). 

In principle, adaptive optics (AO) observations in which images of a PSF are interleaved 
with those of the quasar should be free from geometrical distortion when both are observed on 
the same part of the array. However, our results suggest that the case is not so clear. First, 
the AO observations will suffer from the same nonlinearity issues described above. Second, 
the adaptive correction procedure can be dependent on the object's flux and the details 
of the profile, which can lead to an effective distortion. This underscores the desirability 
of obser ying PSF stars s imult aneously with, and not just close in time to, the quasar; see 



however lAmmons et al.l (120091 ) . 



We conclude that the residual effects of distortion and nonlinearity should be addressed 
by individual near-IR host-hunters for their particular data sets. For the present study, we 
evaluate the residuals based partly on our knowledge of the brightness and proximity of the 
PSF stars. In most cases, we adopt as a criterion that positive residuals are considered 
significant only if they extend beyond a diameter of D > 2.5 FWHM, i.e. a radius r > 
1.25 FWHM, which for the typical frame here means r > 0.55". We explore this further 
with the simulations discussed below. Of course, this particular criterion might not be 
appropriate for data taken under different seeing conditions or with different flux levels. 
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5.4. PSF Fits 

To begin our searcli for tiost galaxies, we modeled each quasar as a point source with 
shape represented by the PSF star images. We determined a two-dimensional (2D) best-fit 
m odel for each quasar using the C program imfitfits provided by Brian McLeod and described 
in Lehar et al. This is the same pr ogram that we used on NICMOS images of low- 

redshift quasars (IMcLeod fc McLeodI 1200 ll ) . We used the 2x magnified images to ensure 
good sampling of the PSF, and extracted an 8" x 8" sub-image for the fitting. Imfitfits 
makes a model by convolving a theoretical point source with the observed PSF, and then 
varying any combination of the parameters defining the background level and the position 
and magnitude of the point source to minimize the sum of the squares of the residuals over 
all the pixels. By subtracting the best-fit model from the quasar image, we can examine the 
result for any residual flux due to an extended component. 

We achieved excellent results with at least one PSF star in most shown in Fig. 

[TUl In a few cases where other sources within the 8" box would bias the fits of the quasar, we 
have simultaneously fit those other sources as either point sources or galaxies. In these cases, 
we subtracted only the fitted quasar component for the figures; the other sources remain for 
comparison. 

Our fitting process necessarily subtracts out any unresolved contribution from the host. 
A logical next step would be to perform a simultaneous fit of a point source plus a model 
galaxy as is commonly done with lower-redshift quasars. Unfortunately, we have found that 
for these data, multicomponent fits give uninterpretable results. The problem is that for 
our data, the likely range of scale lengths for the galaxies are small compared to the seeing 
disk and too little of the galaxy extends beyond the PSF. The result is that running a 
multicomponent fit, whether unconstrained or partially constrained (for example by holding 
fixed the centers, or the centers and the galaxy shape, or the centers and the nuclear flux, 
...), results in the "galaxy" component being turned into a meaningless compact or even 
negative source to improve the fit in the quasar's core, where PSF variations are the biggest. 
One idea to get around this problem is to downweight or mask out the core in the fits. 
However, our tests on real and simulated hosts have shown that the resulting "galaxy" 
is sensitively dependent on the weighting scheme. (We had even seen hints of this with 



the much better resolved hosts in our low-z HST study (iMcLeod &: McLeodI l200ll ).) Thus, 
we have developed a different way to estimate host magnitudes and morphologies via the 
simulations and isophotal diameter analyses described below. 

As another tool for assessing host detection, we have generated one-dimensional (ID) 
radial brightness profiles, measured in circular annuli, of the quasar and PSF images. We 
present the surface brightness profiles in Fig. [TTJ For comparison we also generated a "fully- 
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subtracted" profile by normalizing the PSF to the quasar within the central few pixels and 
subtracting. In most cases, the PSFs are excellent matches to the quasars down to the level 
of the sky noise. 

We caution that the ID profiles need to be interpreted carefully. For example, for our 
ClassicCam image of q0311-5537, the profile alone (see Fig. [11]) looks like those of some host 
detections postulated in the hterature. However, the 2D fit (see Fig. [10]) shows that the 
residual flux is due to PSF mismatch. For the cases where we do have candidate host galaxies, 
we can use the ID proflles to obtain estimates of the host flux. To do this, we subtract a 
fractional PSF proflle that leaves a just-monontonic residual (as any plausible host would 
not decrease in brightness towards its center), and add up the residual light by integrating. 
This technique is necessarily crude, but the data do not warrant more sophisticated fits. 



6. Detection Limits 

In §5] above, we concluded from our PSF-PSF tests that any residual flux from our 
quasar flts outside a diameter of D > 2.5 FWHM is likely signiflcant. In this section, we 
explore this criterion further in two ways: through simulations and through calculations of 
predicted isophotal diameters for galaxies of various types. This second technique is (as far 
as we know) a new and potentially very useful approach for quasar host galaxy studies. 



6.1. Simulated Hosts 

To probe our detection limits we generated suites of fake galaxies and added them to the 
magnifled images of apparently point-like quasars. We selected quasars both brighter and 
fainter than the median for our sample, and images both at, and deeper than, the median 
surface brightness limit. We convolved each model galaxy with the quasar image, added 
the result to the quasar image itself, and reran the analyses with the PSF stars. We also 
duplicated some of the tests with noiseless (Moffat) PSFs having the same FHWM as the 
quasars. Our simulated galaxies included both exponentials (central surface brightness /iq, 
scale length tq) and deVaucouleurs proflles (surface brightness fieff at effective radius r^ff). 
We considered sizes ro,reff of 0.125, 0.25, 0.5, and 0.75", corresponding to 0.88, 1.5, 3.5, and 
5.2kpc at z = 4. These values are similar to the range ob served for z ~ 4 galaxies in the 



Hubble Ultra-Deep Field (HUDF) (lElmegreen et al.l l2007l ) and high-z lensed quasar hosts 



JPeng et al. Ibood l We tested axial ratios 0.2 <h/a<l 



Visual inspection of the residuals supported the validity of our D > 2.5 FWHM criterion; 
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detectable galaxies left residual light outside of that diameter. In terms of flux, we found 
that for the galaxies we tried, the hosts were cleanly visible for an observed Kg — band flux 
ratio F{host) > {nucleus). These hosts leave central (negative) holes in the subtracted 
2D images and have flux in clear excess of the PSF at r 1" in the ID profiles. For 
hosts with F{host) < ^F{nudeus), the detectability by visual inspection depends on the 
size. The hardest galaxies to recover were those whose scale lengths or effective radii were 
< |FWHM, and also the very large r^ff = 0.75" deVaucouleurs galaxies for which too much 
of the galaxy's flux is at low surface brightness. 



6.2. Isophotal Diameter Analysis 

Bolstered by the results from our simulation, we recast our detection criterion from 
Section [5] above as a detection limit in terms of galaxy isophotal diameter Diso, here taken to 
mean the diameter at which the galaxy light drops below the sky noise. In other words, we 
assume that we can detect any hosts that have Disc ^2.5 FWHM for the surface brighteness 
limits of our images. 

To explore the range of galaxies that could be detectable as hosts, we have calculated 



- Disn f or model exponential and deVaucouleur galaxies following the tradition of IWeedman 



( 119861 ) but updated for the currently favored cosmology. We start with exponential and 
deVaucouleurs galaxies covering a range of scale lengths similar to those in §6.1[ We transport 
them to 2; = 4 by applying cosmological surface-brightness dimming and cosmological angular 
diameter distances. We calculate their z = 4 colors and k-corrections by redshifting and 
integrating a spiral galaxy spectral energy distribution template over the filter bandpasses. 
[We have also used bluer and redder templates, but we note that for these data, the k- 
correction is nearly independent of galaxy spectral shape because the observed Kg-hand 
corresponds to rest-frame B.] Finally, we combine these to calculate the isophotal diameters 
for the galaxies given the surface brightness limits of our images. We also compute each 
galaxy's observed magnitude mxs(obs) by integrating the galaxy flux inside the isophotal 
diameter. 

Fig. [T2b shows Diso as a function of observed magnitude m^s (obs) for a range of galaxy 
types and sizes, and for the range of surface brightness limits found for the 2x magnified 
PANIC and NIRI images that we use for host detection. We stress that mxs(obs) represents 
only the fraction of the galaxy's light that falls above the sky noise; it is not simply the 
absolute magnitude adjusted by the cosmological distance modulus. On these plots, only 
galaxies above the 2.5FWHM line are in principle detectable as hosts. One can see that for 
the galaxies considered (i) the faintest visible deVaucouleurs hosts span ~ Imag at a given 
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surface brightness limit; (ii) the faintest visible exponential hosts span ~ 1.5mag; and (iii) 
deVaucouleurs hosts must be relatively brighter to be detected because more of their light 
is hidden in the steeply sloped and unresolved core. 

In Fig. [T2b . we plot these Diso values as a function of L* assuming no evolution in 
the mass-to-light ratio of the stellar population. We adopt for reference a local L* galaxy 
of magnitude M*j^'^^°' = —20.5. Transported to z = 4 such a galaxy would have a total 
magnitude of m*^^^ evolution ~ 23.6, whereas its observed magnitude would be fainter de- 
pending on the surface brightness limit. One can see from the figure how the detectability 
depends on the galaxy scale length. For example, on an image with the median FWHM and 
with the deepest limiting surface brightness (/x = 22.4mag arcsec"^), the intermediate-scale 
(ro = l.Skpc) exponentials are visible at lower luminosity than either the small- or large-scale 
exponentials. The smaller galaxies hide a larger fraction of their flux in the unresolved core. 
The larger galaxies have lower central surface brightnesses at a given luminosity, and their 
relatively more shallow disks do not pop above the surface brightness limit until farther out 
in their profiles. We use these plots to estimate host detection limits for each object. 

We note that there are different ways that one might measure the surface-brightness 
limit of any given image. We have compared our calculated isophotal diameters and ap- 
parent magnitudes from this section with those measured on the images for the simulated 
galaxies discussed above in §6.1[ We find them to be in excellent agreement when the surface 
brightness limit used for the isophotal diameter calculation is that given by the la pixel- 
to-pixel sky noise. This is how we have characterized the surface brightness limits for our 
images in Table [H 



Results 



In a typical image, we are sensitive to field galaxies as faint asm^'^" -23 (m^^ = 24.8), 
with the actual limits dependent upon morphology. To translate this apparent magnitude 
into a corresponding luminosity for a present-day galaxy with the same stellar mass, we need 
to account for luminosity evolution of the stellar population. A reasonable assumption for 
the evolution is that the galaxies undergo ~ 2 mag (i.e a factor of 6) of fading between 
z = A and now, which we infer from the i^-band k-corrections measured for galaxies in 
the HDF-S ( ISaracco et al. II2006I ). This amount of evolution is also expected in the (rest 
frame) visible mas s -to-lig ht ratio according to the stellar population synthesis models of 
Bruzual fc CharlotI (120031 ) for formation redshifts of 2; > 5. If this is the case, our images 
yield galaxies with stellar masses corresponding to a present-day galaxy with luminosity 
< L* in the fields around the quasars. A detailed study of the quasar environments will be 
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presented in Bechtold & McLeod (2009). Here we discuss the results for hosts. 

7.1. Host Limits 

For host galaxies, our detection limits are of course brighter than the limits for galaxies in 
the field. We inspected the radial profiles together with the two-dimensional fits to classify 
detections as y/?/n (likely/maybe/unlikely) with results given in Table HI We looked for 
residuals that extend beyond the sizes of the circles shown in Fig. [TOl and that are not 
likely attributable to nearby (projected on the sky) companions. We further used the ID 
fits to verify that the residuals were plausibly broader than the PSF. We find four likely 
hosts, and note that they are seen on the images that have the best seeing, FWHM < 0'.'4, 
and nearly maximal depth, indicating that we are pushing the limits of detection with these 
images. Other objects might well have hosts lurking just beneath the noise. Three of the four 
likely detections are found on NIRI images. The fuzz associated with quasar q0848 in the 
Ks band was marginally detected in H as well. These four likely hosts are shown in Fig. [T31 

For all of the objects observed with PANIC or NIRI, we estimate the host detection 
limits by applying our Diso ^ 2.5 FWHM criterion using the curves in Fig. [12] and the 
surface brightnesses and FWHMs in Table [H The results are summarized in Table HI where 
we list for each model two possible values for the limit on the host galaxy. 

The "conservative" value represents the most luminous host that would be just visible; 
the galaxies may be luminous but are conspiring to evade detection either by putting too 
much light in their unresolved cores or by having such a large scale length that the middle 
radii are below the sky. The "optimistic" value represents the least luminous host that 
would be just visible. Of course the stellar masses of the galaxies in Table H] could be 
considerably smaller than the straight luminosities indicate. For example, if we allow for 2 
mag of evolution, then the the present-day equivalents would be galaxies lower in luminosity 
by a factor of ~ 6. In that case, a 12L* galaxy in the Table would represent 2L* of stellar 
mass. 

One can see from Table H] that for the depth and resolution of our images, the range 
in upper limits for each type of galaxy typically spans a factor of two. In addition, the 
luminosity limits for deVaucouleurs galaxies are typically double those for exponentials, 
reflecting the fact that the peakier spheroids can hide more light in the unresolved core. In 
general, the least certain limits by the isophotal diameter method will be for those images 
with big FWHM and/or shallow depths, because in these cases galaxy isophotal diameters 
are only weakly dependent on galaxy luminosity-in other words, they fall on the fiat outer 
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parts of the curves in Fig. [T2b . The ClassicCam images were sufficiently insensitive that the 
hmits are not interesting. There are also several PANIC and NIRI images whose conservative 
limits were so large as to be also uninteresting. 



7.2. Host Detections 

For the four likely hosts we can also estimate fluxes from the residuals after fitting. We 
use Fig. [T2b first to identify the kind of galaxy that could yield the observed magnitude and 
isophotal diameter for the surface brightness limit of the image. We then use Fig. [T2b to 
translate that into a possible intrinsic B-band luminosity for the whole galaxy. Note that 
this luminosity is bigger than the luminosity one would calculate simply from applying the 
distance modulus to the observed magnitude; it includes contributions from the inner part 
of the galaxy under the PSF and from the outer part below the sky noise. Finally, we apply 
2 mag of evolution to the luminosity and from it calculate the corresponding Mb that a 
galaxy of the same stellar mass would have today. 

For q0109, the residuals add up to m^'^^'*(obs) ~ 22.7 and they extend to a diameter 
of approximately 2'.'. We use Fig. [T2k (the 22.4mag arcsec"^ depth is appropriate for this 
image) to learn that this galaxy could for example be a large scale-length exponential disk. 
Locating this curve on Fig. [T^ . we find that the same Diso gives a luminosity of ~ SOL* 
with no evolution. Allowing for 2 mag of evolution yields a galaxy with mass corresponding 
to a ~ 5L* galaxy today. For comparison, the object at I'.'S southeast of the quasar has 
m^'^^"(obs) ~ 22.7 and D^so ~ I'-'l- If it is at the redshift of the quasar, it could represent a 
companion with roughly half the mass of the host at a projected separation of about llkpc. 

For the bright residuals in q0234, we measure Diso ~ 2'.'6, while integrating the ID 
residuals gives m]^'^^"(obs) ~ 19.9-bright enough to represent a ~ 60L* exponential (10L*with 
evolution), again with a large scale length. 

A similar analysis for q0848 gives Di^o ~ 1'.'2 and m^'^^"(obs) ~ 19.4 which could be a 
~ 24L* (4L*with evolution) disk. 

For q2047, the residuals extend asymmetrically to the south and possibly represent the 
combined flux of th e host and a companion . This object's classification as a hyperluminous 



infrared galaxy by iRowan- Robinson! (120001 ) (based on sub-mm emission that is likely too 
strong to originate in the quasar's dust torus) would suggest that we are observing a merger. 
Integrating the ID residuals gives m]^^^"(obs) ~ 19.7. The diameter is harder to define 
because of the asymmetry but we adopt 1'.'6 as an estimate, which implies a ~ 40L*-mass, 
intermediate scale length exponential (7L*with evolution). 
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Table m shows that three of the four detections are more luminous than the conservative 
detection thresholds for exponentials and all fall above the optimistic thresholds. We overplot 
the estimates on Fig. [TUfT^ 

If we examine the deVaucouleurs curves in Fig. [12] we find that the residuals for q0109 
and q0234 are too big for their magnitudes observed to be represented by the curves plotted; 
in other words, if they are spheroids, they must have r^ff » 4kpc. On the other hand, the 
residuals for q0848 and q2047 fall on the curves for r^jf = 4kpc and could be spheroids with 
masses of 1-2 times their exponential counterparts. 

For the several quasars with detections listed as "?" the data do not warrant any attempt 
to characterize the magnitudes other than to say that if the hosts are there, they likely lie 
close to the limits listed in Table HI 



7.3. Color of the Host Galaxy of SDSSpJ084811. 52-001418.0 

For one of the quasars with a detected host, SDSSpJ084811. 52-001418.0, we also ob- 
tained a deep H-band image with PANIC. We carried out the PSF estimation and subtraction 
from the nuclear quasar image for the H-band, and found a residual flux, consistent with the 
K-band detection. The color of the host light is H-K=0.8. This color is nominally bluer than 
that of a redshifted spiral galaxy spectral energy distribution, which wo uld have H-K=2. 0. 



The quasar itself has H-K=0.5, as expected for a quasar at this redshift ( IChiu et al.l 120071 ). 



Thus, the host galaxy is redder than the quasar itself, but bluer than a young stellar 
population at z = 4. It could be that the host is experiencing a burst of star-formation, 
as one would expect to occur for a major merger of gas-rich galaxies, and is metal-poor, so 
has a weak 4000Abreak. Another possible explanation, however, is that the host light is 
contaminated by a foreground galaxy. There is no known da mped Ly-a absorber along the 



SDSSpJ084811. 52-001418.0 sight-line flMurphv fc Liskell2004f ) but other intervening absorp- 



tion line systems are no doubt present. Finally, we note that the uncertainty in the H-K color 
is very large, and difficult to estimate. Deeper imaging of a larger sample, or spectroscopy 
of the fuzz, could distinguish among these possibihties. 



7.4. The Local Black-hole/Bulge Relation 

In this section we compare our limits with the local black- hole/bulge relation. While 
the local relation is given for spheroids, we include both exponential and deVaucouleurs 
models in our discussion to allow for the possibility that the stellar mass might be differently 
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distributed at early cosmological times. Wherever necessary we have adjusted values to our 
adopted cosmology, and we have transformed host absolute magnitudes in the literature 
to Mb assuming galaxy rest-frame colors B — \^=0.8, B — R = 1.4, and i? — if = 4, all 
appropriate for a spiral-like stellar population. For an older (elliptical-like) population the 
colors are up to 0.3mag redder, but the uncertainties in these colors are small compared to 
the other uncertainties. 

We have computed rest-frame absolute B magnitude limits for our hosts from the lumi- 
nosities in Table S] assuming 2 magnitudes of evolution. We plot our host limits and black 
hole mass es against the local rela tion in Fig. [TU where t he local popu l ation is represented 



by (i) t he iTremaine et al. I (120021) local galaxies, (ii) the iLauer et al. I (120071 ) local fit, and 



(iii) the lMcLure fc DunlopI (120011 ) local {z < 0.3) luminous quasars, the latter which provide 
high-mass black hole counterparts analogous to those in our sample of high-z quasars. The 
results are similar if we use our OALedd estimates for the black hole mass. We interpret this 
figure as follows. The conservative case occurs where the hosts are all at their maximum 
allowed values, i.e. with limit given by the right-hand bar of each pair and hiding much 
flux in either a very compact core or in the low surface brightness wings of a shallow profile. 
In this case, the exponential and deVaucouleurs distributions are both consistent with the 
local relation for the two magnitudes of evolution assumed; the limits thus do not provide 
interesting constraints on possible evolution in the relation. The optimistic case occurs when 
the left-hand bar in each pair gives the upper limit to the galaxy luminosity. In this case, 
the limits for deVaucouleurs hosts remain uninteresting. On the other hand, hosts with 
exponential profiles would be less massive for a given black hole mass than are the local 
(Tremaine) galaxies, though might yet be consistent with local quasars. However, if the 
evolution correction are more than the two magnitudes assumed, our optimistic exponential 
upper limits would yield hosts lower in mass for a given black hole mass than local luminous 
quasars. The iBruzual fc CharlotI (120031 ) models suggest that evolution in excess of 2 magni- 
tudes between z = 4 and now would be expected if for example the population were younger 
than about 400Myr at z = 4. 

Fig. [15] s hows a similar conapariso n for quasars at various redshifts. The local quasar 
sample is the iMcLure fc Dunlod (I20011) local (z < 0.3) used above. Intermediate-redshift 



Peng et al. (j2006 ) compilation of quasars observed with HST's 



objects {z > 1) are from the 
NIC MOS and having yirial black hole mass est i mates . They include 15 unlensed objects 
from iRidgway et al. I (120011 ) and iKukula et al. I ( l200ll ) and 36 objects (most lensed) from 
their own observations. The points in Fig. (TB] are shown with no evolution correction, which 
explains the (rightwa rd) trend towards more luminous hosts for a given black hole mass 
as redshift increases. iPeng et al. I (120061 ) found that an evolution correction for a passively 
evolving population formed at high redshift pushes the intermediate-redshift hosts leftward 
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beyond the local relation, implying smaller stellar mass relative to the black hole mass 
compared to local galaxies. The constraints from our high-z objects are also dependent on 
the amount of evolution as discussed above. 

In Fig. [161 W6 plot our host limits against rest-frame nuclear 5-band absolute mag- 
nitude, and compa r e them to those in the z < OA quasar host galaxy compilation by 



McLeod fc McLeodI (120011 ). Because the observed Ks directly traces the rest-frame B for 
our quasars, the _B-band absolute magnitude plotted here is independent of the nuclear 
spectral shape and so provides an complementary approach to using black hole mass as a 
tracer of the nuclear engine. We see from the plots that the limits for exponential hosts imply 
galaxies fainter than their low-z counterparts, especially using the optimistic limits (bottom 
bar in each pair) as a bright limit. The same is true for the optimistic deVaucouleurs limits. 
As in the previous discussion, any excess evolution would push the two distributions farther 
apart. An obvious limitation here is that there are few local quasars with luminosities as 
high as those of the z = 4 sample. However, one solid conclusion from Fig. [16] is that 
there does appear to be a maximum allowed host. For the 2 mag of evolution plotted, this 
maximum corresponds to Mb ~ —24 in the conservative limit, or Mb ~ —23 (roughly a 
lOL* galaxy) in the optimistic case. Alternatively, if there are ever independent suggestions 
that the mass limit must be less than that corresponding to an lOL* galaxy, then our results 
would imply that the evolution must be more than the 2 mag assumed. 



7.5. K-band Galaxy Evolution: The K-z relation 

To look at our observations another way, we plot Kg-magnitude versus redshift in Fig- 
ure [T71 The observed i^^-magnitude of a given galaxy varies with redshift because of k- 
corrections, evolution of the galaxy's stellar population, and merging. Here we compare the 
quasar host galaxies of our study with observations of the /^s-magnitude of field galaxies at 
the same redshift. 



We plot observed /To-magnitudes for radio galaxies (ILacy et al.ll2000l : jPe Breuck et al. 



2OO2I : IWillott et al.ll2003l : jPe Breuck et al.ll2006l ). which define the locus of brightest galaxies 
at all redshifts. The locus of radio galaxies is plotted as a solid line given by 

Ks = 17.37 + 4.53/o^ioZ - 0.31{logwzf 



from lWillott et al.l (|2003[ ) . Fainter galaxies f ound as Lyman dropouts (iReddv et al.ll2006ll or 



similar optical selection (jlovino et al.l l2005l : iMcLure et al.l 12006 
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I and 



subsequent spectroscopic redshift measurement are shown as well. We plot Vega magnitudes 
for Kg, and convert from Kab given in the literature by assuming that Kab = Kyega + 1.84. 
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The sharp locus of the radio galaxies is interpreted to indicate a maximum mass for galaxies 
of IO^^Mq, possibly the result of a fragmentation limit in cooling proto-galactic gas clouds. 
Also shown are the expected Kf,-ha.nd evo l ution tracks for elliptical galaxies of various masses 
as c omputed bvlRocca-Volmerange et al.l (120041 ) . The evolution curves for spirals are similar 



see 



Rocca-Volmerange et al.l (j200J)). We see that between z = 0.25 and z = 4, we expect 



a given mass galaxy to undergo approximately 5 magnitudes of evolution at Kg. 

Note that of the four z=4 hosts detected here, three (SDSS010905.8+001617.1, SDSSpJ023446.5^ 
001415.9, and PC2047+0123) are radio quiet. The other, SDSSpJ084811. 52-001418.0, lacks 
definitive radio data. Despite having radio quiet nuclei, two of the three most luminous hosts 
have magnitudes comparable to those of radio galaxies at the same redshift. 

If we interpret the host galaxy detections with the Rocca-Volmerange et al. (2004) 
model for spheroids in the K-z diagram, then we can derive the ratio of black hole mass to 
spheroi d mass: we find a ra tio of 0.016, compared to the local value of 1.4 ± 0.4 x 10~^ seen 
locally (IHaring fc Rixll2004l ). In other words, given the black hole masses we inferred from 
the emission lines, the local relation would imply more massive host ga laxies than we see. 
This is what is predicted by theory (ICrotonll2006l : ISomerville et al.ll2008l ). 



7.6. Malmquist Bias in This Sample 



Lauer et al. I (120071 ) have emphasized that surveys of high redshift quasars suffer from 
such strong Malmquist bias that drawing conclusions about the evolution of host galaxies 
from small samples is problematic. The problem is that we pick targets based on their 
quasar luminosity, then look for the host. Even if there is a correlation between host galaxy 
luminosity and black hole luminosity, the observed host galaxies will be systematically fainter 
than the mean relation derived from spheroid velocity dispersions. This is because for any 
plausible galaxy luminosity function, there are many more faint galaxies than bright ones. 
Thus the intrinsic scatter in the relation causes an excess of low-luminosity galaxies with 
black holes big enough to make the sample. 

To look at this effect quantitatively for our sample, we carried out a Monte-Carlo sim- 
ulation. We chose galaxies randomly from a Schechter galaxy luminosity function, from L/ 
L*= 0.1 to 15, given by 

(l){L)dL = n^{L/UYexp{-L/h*)dL/h* 



with a = —0.46, and L* = 2 x lO^^L© (ISparke fc Gallagheiil2007l ). We assume a galaxy mass 
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to- light ratio M/L=5 (ICappellari et al.ll2006l : iTortora et al.ll2009l ). so that the corresponding 
mass of an L* galaxy is M*^ lO^^Mfr^. Then , we a ssign to each galaxy a black hole with 



mass chosen randomly from the iHaring &: Rixl (120041 ) distribution, in which black holes have 
mass 0.14% +/-0.04% of the galaxy mass. We further assume that to be found as a luminous 
quasar and be eligible for inclusion in our sample, the object must have a black hole mass 
above a certain cutoff, taken to be either logio(MBH/M0) = 8.5 or 9.0. We draw millions of 
galaxies, and count up the fraction of host galaxies as a function of host galaxy luminosity 
for all objects with black hole masses above our threshold. The results are shown in Fig. [TSl 

As expected, we see that for a given value of the black hole cutoff mass, the results will 
be skewed towards lower-luminosity hosts than inferred from the mean value of the black 
hole-bulge relation. If we take a cutoff of logio(MBH/M0) = 8.5, which is at the small end 
of our sample, the effect is modest; most of the hosts would be very close to the expected 
mean value of about 2L*. However, if we restrict our sample to the more massive black holes, 
logio(MBH/M0) = 9, the effect is somewhat larger with the bulk of the contributors in the 
range 4-8L*, skewed from the expected mean value of 7L*. However, even if we include the 
effects of the Malquist bias, we could have detected such hosts for the most massive black 
holes in our sample, those with logio(MBH/M0) > 9.5, if the evolution were at least two 
magnitudes as discussed above. 



8. Discussion 

When we initiated this program, we hoped to test whether or not quasars at z = 4 
followed the same correlation of black-hole mass and host galaxy spheroid mass seen in 
spheroids locally. We knew in principle that if the very luminous nuclei in the high redshift 
objects were hosted by proportionately massive spheroids, they would be relatively straight- 
forward to detect at K, with the best ground-based seeing. We have taken a very conservative 
approach to reducing the data, to looking for host galaxy detections, and to estimating host 
galaxy fluxes and limits. 

We explored more than one way to compare our observed detection of hosts and limits 
on host galaxies to the low redshift spheroids. This comparison is complicated for a number 
of reasons, primarily the fact that we do not measure spheroid velocity dispersion or mass 
directly, but must infer the host galaxy properties from the emitted starlight. Nonetheless, 
our data indicate that the host galaxies of some z=4 quasars in our sample are fainter than 
we expect from the low-redshift correlations. 

We note that for the mean local values MsH/Mgai = 0.14% and M/L=5, we expect 
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black holes with logio(MBH/M0) = 9.5, 10, and 10.5, to have host galaxies with Mgai = 
2,7, and 22xlO^^M0, or L/L* = 20,70, and 220 respectively. However, such galaxies are 
implausibly large, and do not have present-day analogues. For example, as shown in Fig. 
[T71 the upper mass envelope for radio galaxies corresponds to lO^^M©, or L ^ lOL*. This 
suggests that the MsH/Mgai correlation must plateau at Mgai ~ lO^^M©. 



9. Summary 

We observed 34 high redshift quasars in the near-IR to search for their host galaxies. 
Our conclusions are the following. 

1. We found that to characterize the PSF and subtract the nuclear quasar light properly, 
we had to account for geometric distortions in the camera optics, and non-linearity in the 
detectors, beyond what is normally corrected for in standard pipeline reductions. We caution 
that small uncertainties in the linearity correction, and the use of PSF stars which are 
significantly brighter than the quasar, can lead to undersubtraction of the nuclear PSF, and 
spurious host galaxy detections. 

2. We derived black hole masses for the quasars in the sample from the profile of the C 
IV emission line, but noted several cases where the profile appears to include a narrow and 
broad component. Low signal-to- noise spectra could easily miss the broad component, with 
the result that the black hole mass is underestimated. The black hole masses range from 
lO*'^ to IO^'^'^Mq. These quasars are very luminous, and so rare a s to b e not represented 



i n som e models for quasar evolution (e.g. iKauffmann &: Haehneltl (120001 ): iDi Matteo et al. 



( I2OO8I )). They are more 



uminous than the k nee in the quasar luminosity function, and are 



probably peak emitters (iHopkins et al.ll2008l ) that are undergoing a major merger 



3. Accretion rates were derived from the observed K photometry, which directly samples the 
rest-frame B for the sample quasars. The median accretion rate of the sample corresponds 
to an Eddington fraction of Lhoi/LEdd = 0.41 ± 0.3, consistent with the findings for other 
samples such as the Sloan quasars, and for low redshift quasars whose host galaxies have 
been studied by a number of authors. 

4. We estimated the K-band magnitudes of the host galaxies of the quasars in our sample with 
a new method which takes into account the isophotal diameter of the galaxies as a function of 
redshift, as well as the surface brightness limit of the images. We explored parameter space 
by considering galaxies with exponential and deVaucouleur radial light distributions, and 
quantified the dependance of derived host galaxy properties on assumed galaxy properties. 
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5. We detected host galaxies for 4 quasars, at least three of which are radio quiet. T/ie 
detections all occurred on our deepest, sharpest images. The K-band luminosities of the 
hosts are consistent with massive galaxies at the redshifts of the quasars. For one object 
with H-band data, the K-H color is bluer than expected for a star-forming galaxy at the 
quasar redshift, but the uncertainties in the color are large. 

6. We interpreted the detected hosts and limits on host luminosity in several ways, taking into 
account expected evolution of the stellar populations. We find that if the hosts are already 
spheroids at early times, then their black- hole/bulge relation could well be consistent with 
that for local galaxies and luminous quasars; our limits arc weak in the case of very compact 
or very extended spheroids. On the other hand, if the hosts are exponential disks, they 
likely have less stellar mass for a given black hole mass than would be inferred from the 
extrapolation of the local relation to high black hole masses, but they could contain as much 
stellar mass as the spheroidal component of local luminous quasars. Any S-band rest frame 
luminosity evolution in excess of the 2 magnitudes assumed would make any evolution in the 
black-holc/bulge relation stronger; such would be the case for a stellar population younger 
than -.'400Myr. 

7. If we interpret the K-magnitudes of the detected hosts with models for the evolution of 

massive spheroids we find that the ratio of black hole mass to spheroid mass for the 4 detected 
hosts is approximately 0.02, compared to 1.4 ±0.4 x 10^^ observed in local spheroids. Several 
authors have pointed out that the Malmquist bias inherent in any study that looks for hosts 
in very luminous, rare quasars will overestimate the black hole mass to spheroid ratio, given 
the likely scatter in the correlations, and the fact that faint galaxies outnumber bright ones. 
We made a rough estimate of the Malmquist bias in our sample through a Monte Carlo 
simulation. We find that our detection rate is inconsistent with what we should have seen, 
had the z=4 quasars followed the same relation between black-hole mass and spheroid mass. 
Instead, the host galaxies in the past appear fainter (and by assumption less massive) than 
host galaxies today. This conclusion depends on uncertain assumptions for the scatter in the 
black-hole mass-spheroid correlation, the mass-to-light ratios of galaxies, and the evolution 
of the spectral energy distribution of spheroids. However, our results are in broad agreement 
with semi-analytical models for the growth of black holes and merger-induced activity. 
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Fig. 1. — Magnitude and redshift distributi on for the 34 quasars observed here (squares) 
and quasars from the Sloan Quasar Catalog (jSchneider et al. 1120071 ). 
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Fig. 2. — New Keck DEIMOS spectra of five quasars in the region around CIV in the rest- 
frame used for estimating black hole masses. The spectrum for [VH] 2 125-4529 also provided 
an improved redshift. 
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Fig. 3. — Black hole masses as a function of z for the 34 quasars observed here (filled 
squares) shown on top of those determined for ^ 1600 SDSS quasars in this redshift range by 



Shen et al.l (120081 ) . For several objects the CIV line profile shows a distinct broad component 



whose FWHM yields black hole masses (open squares) an order of magnitude larger than 
the overall FWHM would indicate. 
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Fig. 4. — Sample PANIC image showing a, z — 4.0 quasar "q" and PSF stars "a" and "c." 
The surface brightness hmit in the center is K — 22mag arcsec"^, which is shghtly better 
than average, but at the seeing is shghtly worse than average. At this redshift, the field 
shown is ~ IMpc across. 
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Fig. 5. — Effects of distortion corrections and recentering for star "c" from Fig. HI Top 
left: frames have been aligned to the quasar centroid before combining, and no distortion 
correction has been performed. The star "c" image is stretched. Top right: the nominal 
distortion correction has been applied before aligning to the quasar centroid and combining, 
which has tightened up the star "c" image; however, it is still broader than the image of 
the quasar taken from the same frame. Bottom left: the distortion-corrected frames have 
instead been aligned to the centroid of star "c" before combining, effecting a second-order 
distortion correction for that star. The star "c" image now has the same FWHM and shape 
as the quasar did in the image aligned on the quasar centroids. 
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Fig. 6. — Improvement of fits for the quasar shown in Fig. H] and [5] after creation of "re- 
centered" PSF frames for stars "a" and "c." The better fit with star "a" is hkely due to its 
relative proximity to the center of the frame. 
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Fig. 7. — Nonlinearity data for PANIC camera based on exposures of internal calibration 
lamps with three possible prescriptions for linearity corrections. Exposure times are generally 
chosen to keep quasars (and PSF stars) below about 15,000 ADU. At this level the possible 
corrections vary by < 0.5%. 
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Fig. 8. — Effects of imperfect linearity corrections on PSF fits. Top row: simulated stars 
with the same FWHM and noise characteristics as our observed quasars and having relative 
flux 1, 10, and 40 times. Nonlinearity has intentionally been applied according to the middle 
curve from Fig. Bottom three rows: results of correcting and subtracting these stars 
from each other via different schemes. In rows 2 (PSFs normalized to the central pixel 
before subtraction) and 3 (normalization based on a fit that minimizes residuals) we have 
performed no linearity correction. Residuals are obvious and the details depend on the 
normalization. In row 4, we have applied a linearity correction that differs from the one used 
to generate the stars by < 0.5%. The spurious residuals are only obvious with the brighter 
PSFs, and do not extend beyond a diameter of 2.5*FWHM. 
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Fig. 9. — Radial profiles for the simulated stars described in §5.21 and shown in Fig. [HI The 
stars that have not been corrected for linearity. For p40, whose flux puts it into the nonlinear 
regime, the profile is obviously different. For plO and pi, the more subtle nonlinearity falsely 
suggests a host contribution in the profile, though the residuals are unphysically compact. 
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Fig. 10. — Closeups of quasars (left image in each pair) and residuals after PSF fits, ~ 8" 
on a side for targets observed with NIRI. Images are shown with North up, East to the left. 
The circles have diameter D = 2.5FWHM. The obviously bad fit for q0241H resulted from 
telescope mirror support problems, and illustrates the effects of distortion in an extreme 
case. The q0120 image shows the much more subtle but typical effect of residual distortion 
and nonlinearity. In the case of q0234, significant residuals are obvious. 



-45- 




Figure [TU] cont'd-targets observed with PANIC. In most cases the fit to the core is excellent. 
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Figure fTOl cont'd-targets observed with PANIC 
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Figure [To] cont'd-t ar gets observed with ClassicCam. Fits are generally poorer than those 
for NIRI and PANIC. As described in the text, for the small ClassicCam field of view, PSF 
stars were not visible on the same frames as the quasar and had to be obtained in separate 

observations interleaved with the quasar exposures. In addition, the ClassicCam images 
generally have broader PSFs and shallower depths. 
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Fig. 11. — Radial profiles labeled with instrument names (Cl,C2=ClassicCam, N=NIRI, 
P=PANIC). Quasar data points are plotted with standard-of-the-mean errors. PSFs are 
shown as solid lines normalized to the central pixel. Dotted lines show positive parts of 
residuals after subtraction. Vertical bars mark the radius corresponding to a diameter as 
D = 2.5FWHM discussed in the text. 
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Fig. [n] cont'd. 
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Fig. 12. — Isophotal diameter predictions v. observed magnitude rri^^°'{ohs) (a) and lumi- 
nosity L/L* (b) for model galaxies at z = 4. The observed magnitudes include only the 
galaxy flux that is inside the isophotal diameter. The luminosity plots assume no evolution; 
with a more realistic 2 magnitudes of evolution the luminosity axis labels will decrease by 
a factor of ~ 6. The quantity r represents ro and Vgff for exponential and deVaucouleurs 
laws respectively. The three vertical panels in each set represent surface brightness limits 
appropriate for the range of our NIRI/PANIC observations. Dashed horizontal lines show 
the median FWHM and 2.5 FWHM. As described in §6.2^ we set detection limits for the 
various host galaxy types according to the criterion Diso ^ 2.5 FWHM. 
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Fig. 13. — Closeups of likely hosts. Boxes are 4"on a side with North up, East left, and 
circles as in Fig. [TDi Top row: NIRI targets. Bottom row: PANIC object q0848 in K and 
H. 
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Fig. 14. — Comparison with the l ocal black-hole/ 



resented by t he op en circles from iTremaine et al. 



j ulge r elation. Local galaxies are rep- 
plus the d i agona l line fit from 



Lauer et al. I (120071 ). Local luminous quasars from 



McLure fc DunlopI (120011 ) are shown as 



filled squares. For each of our objects, the two connected vertical bars mark the optimistic 
(left end) and conservative (right end) upper limits on luminosity for the range of galaxy 
scale lengths represented in Fig. [121 Limits are plotted for (a) exponential hosts and (b) 
deVaucouleurs hosts. Triangles indicate the host luminosity estimates for the 4 host galax- 
ies detected. Our host magnitudes and limits have been corrected for 2 mag of evolution. 
Data from ClassicCam are not plotted, for clarity. All values from the literature have been 
adjusted to Hq = 70 km s~^ Mpc~^. 
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Fig. 15. — Comparison with quasars at different redshifts. Local luminous quasars (filled 
squ ares) and the diago nal line are the same as in Fig. [TH Intermediate redshift quasars from 
the iPeng et al. I ( 120061 ) compilation are shown as circles. For each of our objects, limits are 
plotted as in Fig. [T3] except that here we have made no correction for evolution. Data from 
ClassicCam are not plotted, for clarity. 
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Fig. 16. — Nuclear abs olute magnitude v. galaxy magnitude. Squares show local quasars 
from the compilation of iMcLeod fc McLeodI (120011 ) adjusted to Hq = 70 km s~^ Mpc~^. For 
each of our high-z quasars, the two vertical bars mark the optimistic (bottom end) and 
conservative (top end) absolute magnitude limits for the range of galaxy scale lengths rep- 
resented in Fig. [12] and assuming 2 magnitudes of evolution. Triangles indicate the four 
detected hosts at z=4. For clarity, the data from ClassicCam are not plotted. 
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2 4 6 

Redshift 

Fig. 17. — K-z Diagram for Radio Galaxies, Field Galaxies, and z — 4 Quasar Hosts. The 
filled triangles are the four detected host galaxies of this survey; open triangles are optimistic 
limits listed in Table 4. Open circles are radio galaxies (Lacy et al. 2000, De Breuck et al. 
2002, Willott et al 2003, De Breuck et al 2006). X's are optically-selected galaxies (Reddy 
et al. 2006, McLure et al 2006), and dots are galaxies in the VIRMOS deep galaxy survey 
(lovino et al. 2005, Temporin et al. 2008). The sohd hne is the fit to the locus of radio 
galaxies, as given by Willott et al. (2003). The dotted, dashed, long-dashed and dot-dashed 
lines are evolutionary tracks for elliptical galaxies of mass IO^^Mq, 10^^Mq,10^-'^Mq, and 
1O^°M0, respectively, as computed by Rocca-Volmerange et al. (2004). All magnitudes are 
Ks in the Vega system. The detected quasar hosts have Kg magnitudes consistent with those 
of massive galaxies at their redshifts. 
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Fig. 18. — Results of Monte Carlo simulation of Malmquist bias for M/L=5. Curves represent 
the fraction of bright galaxies at each luminosity expected to have black hole masses above 
the cutoff value shown when drawn at random from a Schechter function and assigned a black 
hole mass according to the present-day Magorrian relation and its scatter. The vertical bars 
show the corresponding galaxy mass predicted from the Magorrian relation with no scatter. 
For the bigger cutoff black hole mass, the Malmquist bias is more severe, in that there is a 
larger contribution from galaxies with lower-luminosities than that inferred from Magorrian 
relation mean. 
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Table 1. Sample and Observations 



Time FWHM SBlim'^ 



Name 




z 


TAB 


rcP 


Date 


Inst'' 


Filter 


(min) 




{") 


magi 


SDSSpJ010905.81+001617.1 


3 


680 


20 


50 


FanOl.SS 


2003 


Sep 11 


N 


Ks 


250 





35 


22 


62 


SDSSJ012019.99+000735.5 


4 


101 


19 


96 


SS 


2003 


Nov 05 


N 


Ks 


125 





32 


22 


41 


SDSSpJ012650.77+011611.8 


3 


660 


19 


20 


Fan99 


2003 


Nov 03 


N 


Ks 


125 





41 


22 


36 


SDSSpJ023446. 58-001415. 9 


3 


600 


20 


17 


FanOl.SS 


2003 


Sep 12 


N 


Ks 


245 





39 


22 


64 


BRI0241-0146 


4 


053 


18 


50 


SL96,ACS 


2005 


Jan 26-27 


P 


Ks 


143 





55 


21 


67 


BRI0241-0146 


4 


053 


18 


50 


SL96,ACS 


2003 


Sep 10 


N 


H 


219 





45 


23 


28 


BRJ0301-5537 


4 


133 


19 


20 


POl 


2002 


Sep 22 


CI 


Ks 


108 





82 


20 


84 


BRJ0311-1722 


4 


039 


18 


30 


POl 


2002 


Sep 21 


CI 


Ks 


129 





38 


20 


83 


SDSSpJ031427.92+002339.4 


3 


680 


20 


55 


FanOl.SS 


2003 


Nov 06 


N 


Ks 


125 





57 


22 


30 


SDSSpJ032226. 10-055824.7 


3 


957 


20 


07 


SS 


2004 


Jan 05 


P 


Ks 


135 





41 


21 


76 


SDSSpJ034109. 35-064805. 1 


4 


070 


20 


48 


SS 


2004 


Jan 04 


P 


Ks 


135 





41 


21 


95 


BR0401-1711 


4 


227 


19 


60 


POl 


2005 


Jan 28 


P 


Ks 


94 





43 


21 


58 


SDSSpJ040550.26+005931.2 


4 


050 


20 


22 


FanOl 


2004 


Jan 02 


P 


Ks 


128 





35 


21 


81 


BRJ0529-3552 


4 


172 


18 


80 


POl 


2003 


Jan 20 


C2 


Ks 


108 





42 


20 


45 


PMNJ0741-5305 


3 


743 


21 


90 


LOl 


2004 


Jan 06 


P 


Ks 


135 





47 


21 


87 


SDSSpJ084455. 08+001848.5 


3 


690 


19 


90 


S01,SS 


2004 


Jan 05 


P 


Ks 


142 





45 


21 


78 


SDSSpJ084811. 52-001418.0 


4 


124 


19 


24 


SS 


2004 


Jan 02 


P 


Ks 


162 





36 


21 


98 


SDSSpJ084811. 52-001418.0 


4 


124 


19 


24 


SS 


2004 


Jan 04 


P 


He 


199 





41 


22 


82 


SDSSpJ085430. 18+004213.6 


4 


084 


20 


29 


SS 


2004 


Jan 03 


P 


Ks 


155 





50 


21 


70 


SDSSJ094822.97+005554.3 


3 


878 


20 


30 


SS 


2004 


Mar 09 


P 


Ks 


196 





39 


22 


00 


SDSSJ094822.97+005554.3 


3 


878 


20 


30 


SS 


2005 


Jan 28 


P 


He 


203 





70 


22 


81 


SDSSJ095755.63-002027.5 


3 


753 


19 


36 


SS 


2005 


Jan 26 


P 


Ks 


250 





46 


22 


16 


SDSSJ095755.63-002027.5 


3 


753 


19 


36 


SS 


2005 


Jan 27 


P 


He 


202 





42 


22 


32 


SDSSJ100151. 58-001627.0 


3 


674 


19 


31 


SS 


2004 


Mar 11 


P 


Ks 


169 





57 


21 


85 


SDSSpJ101832. 46+001436.4 


3 


826 


20 


22 


SS 


2004 


Mar 10 


P 


Ks 


209 





40 


21 


87 


SDSSpJ102043. 82+000105.8 


4 


207 


20 


53 


SS 


2004 


Mar 08 


P 


Ks 


182 





43 


21 


88 


SDSSJ135134.46-003652.2 


4 


010 


20 


13 


SS 


2004 


Mar 11 


P 


Ks 


209 





51 


21 


92 


SDSSJ151618.44-000544.1 


3 


735 


20 


23 


SS 


2004 


Mar 08 


P 


Ks 


148 





41 


21 


86 


SDSSpJ152443. 19+011358.9 


4 


095 


20 


14 


SDSS2 


2003 


Jul 17-18 


P 


Ks 


213 





58 


21 


67 


SDSSpJ154014.57+001854.7 


3 


830 


20 


41 


S01,SS 


2004 


Mar 10 


P 


Ks 


169 





50 


21 


89 


SDSSpJ161926. 87-011825. 2 


3 


840 


20 


17 


FanOO 


2004 


Mar 09 


P 


Ks 


155 





44 


21 


90 


SDSSpJ165527.61-000619.2 


3 


990 


20 


05 


FanOO 


2003 


Jul 19 


P 


Ks 


157 





64 


21 


67 


[HB89] 2000-330 


3 


773 


17 


6 


HB,WFPC2 


2002 


Sep 22 


CI 


Ks 


117 





82 


20 


72 


PC2047+0123 


3 


799 


19 


46 


S91 


2003 


Sep 10-11 


N 


Ks 


154 





35 


22 


46 


[VH95]2125-4529 


3 


763<^ 


19 


50 


HV96 


2003 


Jul 18 


P 


Ks 


146 





56 


21 


54 


[VCV96]Q2133-4625 


4 


180 


21 


23 


VV,HV96 


2003 


Jul 17 


P 


Ks 


243 





66 


21 


84 


BR2212-1626 


3 


990 


18 


80 


SL96 


2002 


Sep 21 


CI 


Ks 


117 





42 


20 


69 


BR2212-1626 


3 


990 


18 


80 


SL96 


2003 


Jul 19 


P 


He 


157 





57 


21 


89 


BRJ2317-4345 


3 


943 


18 


80 


POl 


2002 


Sep 20 


CI 


Ks 


117 





46 


20 


69 



^"S ources for z^tau - SS ( SPSS D R6); Fan99.00.01 llFan et al.lll999 l. l2000l.l200ll): LOl llLandt et al. llioOlh: POl llPeroux et an200lh : 
SOI llSchneider et al.ll2"oOlb ; SL96 lIStorrie-Lombardi et al.lll996l); HB llHewitt fc Burbidgelll989ll : VV llVeron-Cettv fc Veronill996l . 
I2OOII) : ACS,WFPC2 archival HST images measured here; HV96 iHawkins fc Veronlll99i rand r^B = R + 0.3. 

''Cl,C2=ClassicCam with scale O'.'115/pix, O'.'95/pix; N=NIRI with scale O'.'116/pix; P=PANIC with scale 0^'125/pix 



'^Surface brightness limit of final image (Vega mags) derived from la pixel-to-pixel variation. Effective limits for the magnified 
images used with PSF fits are ~ 0.4mag brighter. 



Improved redshift from Keck spectrum in this paper 



Table 2. Radio Properties 

















RL or 


Name 


z 


E(B-V) 




i rcf*^ 


(mjy) Flux rof^ 


t 


R (2o-)<= RQ? 



SDSSpJ010905. 81+001617.1 


3, 


.680 


0, 


,026 


20, 


,31 


SDSSJ012019.99+000735.5 


4, 


,101 


0, 


,037 


19, 


,82 


SDSSpJ012650.77+011611.8 


3, 


,660 


0, 


,026 


19, 


,17 


SDSSpJ023446. 58-001415. 9 


3, 


,600 


0, 


,022 


19, 


,88 


BRI0241-0146 


4, 


,053 


0, 


,028 


18, 


,24 


BRJ0301-5537 


4, 


,133 


0, 


,010 


18, 


,97 


BRJ0311-1722 


4, 


,039 


0, 


,032 


18, 


,01 


SDSSpJ031427.92+002339.4 


3, 


,680 


0, 


,100 


20, 


,04 


SDSSpJ032226. 10-055824.7 


3, 


,957 


0, 


,053 


19, 


,86 


SDSSpJ034109. 35-064805. 1 


4, 


,070 


0, 


,053 


20, 


,09 


BR0401-1711 


4, 


,227 


0, 


,026 


19, 


,33 


SDSSpJ040550.26+005931.2 


4, 


,050 


0, 


,444 


19, 


,08 


BRJ0529-3552 


4, 


,172 


0, 


,031 


18, 


,51 


PMNJ0741-5305 


3, 


,743 


0, 


,227 


21, 


,08 


SDSSpJ084455. 08+001848.5 


3, 


,690 


0, 


,035 


19, 


,74 


SDSSpJ084811. 52-001418.0 


4, 


,124 


0, 


,030 


18, 


,93 


SDSSpJ085430. 18+004213.6 


4, 


,084 


0, 


,049 


19, 


,84 


SDSSJ094822.97+005554.3 


3, 


,878 


0, 


109 


19, 


,85 


SDSSJ095755.63-002027.5 


3, 


,753 


0, 


,033 


19, 


,28 


SDSSJ100151. 58-001627.0 


3, 


,674 


0, 


,037 


19, 


,18 


SDSSpJ101832. 46+001436.4 


3, 


,826 


0, 


,043 


19, 


,98 


SDSSpJ102043. 82+000105.8 


4, 


,207 


0, 


,040 


19, 


,75 


SDSSJ135134.46-003652.2 


4, 


,010 


0, 


,034 


19, 


,79 


SDSSJ151618.44-000544.1 


3, 


,735 


0, 


,057 


19, 


,89 


SDSSpJ152443. 19+011358.9 


4, 


,095 


0, 


,055 


19, 


,88 


SDSSpJ154014.57+001854.7 


3, 


,830 


0, 


,092 


20, 


,06 


SDSSpJ161926. 87-011825. 2 


3, 


,840 


0, 


117 


19, 


,69 


SDSSpJ165527.61-000619.2 


3, 


,990 


0, 


,298 


19, 


,54 


[HB89] 2000-330 


3, 


,773 


0, 


130 


17, 


,04 


PC2047+0123 


3, 


,799 


0, 


105 


19, 


,07 


[VH95]2125-4529 


3, 


,763 


0, 


,022 


19, 


,15 


[VCV96]Q2133-4625 


4, 


,180 


0, 


,025 


21, 


,04 


BR2212-1626 


3, 


,990 


0, 


,028 


18, 


,49 


BRJ2317-4345 


3, 


,943 


0, 


Oil 


18, 


,57 



SS 


< 0.110 


FIRST 


18.80 


0.90 


Q 


SS 


< 0.112 


FIRST 


18.78 


0.72 


Q 


Fan99 


0.210 


COl 


18.09 


0.43 


Q 


SS 


< 0.108 


FIRST 


18.82 


0.73 


Q 


ACS 


< 0.150 


FIRST 


18.46 


0.13 


Q 


POl 




[NED] 






7 


POl 


< 2.500 


NVSS 


15.41 


1.34 


7 


SS 


< 0.112 


FIRST 


18.78 


0.81 


Q 


SS 


< 2.500 


NVSS 


15.41 


2.08 


7 


SS 


< 2.500 


NVSS 


15.41 


2.17 


7 


POl 


< 2.500 


NVSS 


15.41 


1.87 


7 


FanOl 


< 2.500 


NVSS 


15.41 


1.77 


7 


POl 


< 2.500 


NVSS 


15.41 


1.54 


7 


LOl 


44.000 


G93 


12.29 


3.51 


L 


SS 


5.830 


NVSS 


14.49 


2.10 


L 


SS 


< 2.500 


NVSS 


15.41 


1.71 


7 


SS 


< 2.500 


NVSS 


15.41 


2.07 


7 


SS 


2.320 


FIRST 


15.49 


1.75 


L 


SS 


< 0.149 


FIRST 


18.47 


0.63 


Q 


SS 


< 0.146 


FIRST 


18.49 


0.58 


Q 


SS 


< 0.151 


FIRST 


18.45 


0.91 


Q 


SS 


2.000 


FIRSTc 


15.65 


1.64 


L 


SS 


< 0.148 


FIRST 


18.47 


0.83 


Q 


SS 


< 0.045 


COl 


19.77 


0.35 


Q 


SS 


< 0.143 


FIRST 


18.51 


0.85 


Q 


SS 


< 0.151 


FIRST 


18.45 


0.94 


Q 


FanOO 


< 0.025 


COl 


20.41 


0.01 


Q 


FanOO 


0.010 


COl 


21.40 


-0.74 


Q 


WFPC2 


570.00 


W90 


9.51 


3.01 


L 


S91 


0.105 


S92 


18.85 


0.09 


Q 


HV96 




[NED] 






7 


HV96 




[NED] 






7 


SL96 


< 2.500 


NVSS 


15.41 


1.54 


7 


POl 




[NED] 






7 



''Corrected for reddening 

^S (SDSS DR6); Fan99.00 .01 llFan et al.ll999l.l2000ll200ll) : LOl l lLandt et al. I2OOII) computed from : POl jPeroux et all 

12001*) computed from R; S91 fSchncider et al. "l99l') computed from Agl450: SL96 l lStorrie- Lombard! et al.lll996l ') computed 
from AB7000A; HV96 (Hawkins fc Veroniil996i) : ACS.WFPC archival image measured here 

■^Detections or Icr upper limits 

■^FIRST and NVSS were accessed through NED; COl l lCarilli et al.l I2OOII') : G93 llCriffith fc Wriehtl Il993h : S92 
dSchneider et al. \\l99i) : W90 llWright fc OtrupcekllTgOol) 

"Radio-to-optical flux ratio calculated from detections, or from 2(t limits for nondetections 
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Table 3. Black Hole Properties 
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SL96 


2739/6185^ 


C02 


-28.5 
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POl 


-27.4 
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''Reddening-corrected AB magnitude at 1450A(1 + 2) 

^Sources for spectra and p hotome try: i (d erived from i^^ in Table[2ll: C02 llConstantin et al.l2002h: D03 dPietrich et al.'2003'): Fan99,00.01 
llFan et al.l[l9 99. 200oi. l200lh : KOI llKuhn et al. 20o 3); LO l l|Landt et al. Il200lh : P Ol jPeroux et al.ll200l|): SOI. 92 ( Schneider et al. 200ll : 
ISchneider etal. 111993 '): SL96 lIStorrie-Lombardi et anfl996h : SS (SDSS DR6); SWOO lIStorrie-Lombardi fc 'VVoif3l2000l) : Keck this paper 

'^FWHM measured by us from spectra. For objects with two values Usted we measured both a broad and a narrow component. 

'^Quasar absolute magnitude Mg (Vega system) computed from the observed K-band magnitudes in Table |3] 

''Black hole mass logioCM^n/TsAQ) (derived from FHWMcjv and AB1450) and Eddington fraction Li,oi/ L^dd (from mass and Mg). 
''o9io(Mbh/Mq) calculated from quasar magnitude Mg assuming L/L^^^ = 0.4 

^Notes: q0405 and ql655 profiles have broad wings; q0529 spectrum is noisy; q2212 is gravitationally lensed llWarren et al.ll200ll) 
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Table 4. Results from Kg Imaging 
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'^Range of upper limits for exponential liosts with no evolution-2mag evolution would lower L's by 6x 
''Range of upper limits for dcVaucouleurs models with no evolution— 2mag evolution would lower L's by 6x 
•^Observed K magnitude of galaxy corresponding to min Lea;p/L* in column 4 
■^Estimate of host luminosity for disks with and without 2mag of evolution 
''Host absolute magnitude in the rest frame assuming 2 mag of evolution 



